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Experimental determination of the state- dependent enhancement of the 
electron-positron momentum density in solids 

J. Laverock, T. D. Haynes, M. A. Alam and S. B. Dugdale 

H. H. Wills Physics Laboratory, University of Bristol, 

Tyndall Avenue, Bristol BS8 ITL, United Kingdom 

The state-dependence of the enhancement of the electron-positron momentum density is inves- 
tigated for some transition and simple metals (Cr, V, Ag and Al). Quantitative comparison with 
linearized muffin-tin orbital calculations of the corresponding quantity in the first Brillouin zone 
is shown to yield a measurement of the enhancement of the s, p and d states, independent of 
any parameterizations in terms of the unscreened electron density local to the positron. An em- 
pirical correction that can be applied to a first-principles state-dependent model is proposed that 
reproduces the measured state-dependence very well, yielding a general model for the enhancement 
of the momentum distribution of positron annihilation measurements, including those of angular 
correlation and coincidence Doppler broadening techniques. 



I. INTRODUCTION 

Positron annihilation is one of the key tools in mod- 
ern investigations of the Fermi surface (FS) of solids [1], 
alongside quantum oscillatory techniques, Compton scat- 
tering and angle-resolved photoemission. Unlike other 
methods, however, the positron probe itself plays a cru- 
cial role in the measured distribution, preferentially anni- 
hilating with those electrons that are most able to screen 
its charge. In a typical metal, free from vacancy-type de- 
fects, the electrons that are most readily able to screen 
are, of course, those at the FS, advantageously leading to 
an enhancement of the signal contributed from electrons 
at the FS. Attempts to account for this enhanced contri- 
bution have, for the most part, relied on detailed stud- 
ies of the electron-positron interaction within the jellium 
model [2] |3] , which is now essentially well- understood [4] . 
However, such schemes are yet to achieve good agreement 
with experiment when applied to a wide range of metallic 
systems. Here, we consider this problem from an exper- 
imental perspective, measuring the state-dependent en- 
hancement factor for some simple elemental metals, and 
present a phenomenological (and empirical) correction to 
the work of Barbiellini, Alatalo and their co-workers [5l|6] 
that offers excellent agreement with experiment. 

The complex many-body interaction between the 
positron and the electron gas has been intensively studied 
for many years [7 . When the positron enters a homoge- 
neous electron gas, the attractive Coulomb interaction 
polarizes the electron gas in the vicinity of the positron, 
leading to a cusp in the unscreened electron density at 
the positron's position and the associated enhancement 
of the partial annihilation rate of those electrons that 
screen the charge. The theory of Kahana [2] predicted 
a momentum-dependent enhancement, in which the en- 
hancement increases towards the Fermi momentum, /cp, 
and corresponds to the increased capability of electrons 
near the Fermi level to screen the positron's charge, com- 
pared with lower-lying electron states. However, the in- 
homogeneity of the electron gas in real lattices can have 
a strong influence on the enhancement, even hiding the 



Kahana-like momentum dependence [5j [7] . 

It is worth pointing out that when considering en- 
hancement there are actually two separate, but related 
issues. Firstly, in the context of calculating the correct 
positron lifetimes in solids, the enhancement of the to- 
tal electron density needs to be properly described in 
order to calculate the positron annihilation rate. Sec- 
ondly, a description of the enhancement is needed when 
calculating the two-photon momentum densities (which 
are the focus of the current paper). The former problem 
is easier because the contact density can be parameter- 
ized in terms of the local electron and positron densi- 
ties (using the many-body results for jellium), but the 
latter is a more difficult problem since in the frame- 
work of density functional theory there is no formally 
exact way to calculate the two-photon momentum den- 
sity [HI E] (and as such all models in the literature are, 
in practice, empirical). Local density parameterizations, 
in which the enhancement is parameterized as a function 
of the unscreened electron density, n, at the positron, 
have been introduced to account for the inhomogeneity 
of real systems. In these approaches, the enhancement is 
usually expressed in terms of the electron gas parameter, 
Ts = (3/47rn)^/^, of which it is a monotonically increas- 
ing function for typical crystallographic electron densi- 
ties. Some popular choices are the expressions of Arpo- 
nen and Pajanne [3 , based on boson formalism and pa- 
rameterized by Barbiellini and co-workers [10 , and those 
of Boronski and Nieminen (BN) which are based upon an 
interpolation of Fermi liquid results due to Lantto [TT]. 
Jarlborg and Singh (JS) have used a local-density ap- 
proach to solve a two-body electron-positron Schrodinger 
equation inside a spherical correlation cell that yields 
good agreement with transition metals and their alloys 
for both momentum densities [T2^ and positron lifetimes 
[13 , and is a common choice to describe the enhance- 
ment of the momentum distribution in metals [9l IT4] . 
More general parameterizations have been proposed [15]- 
18 that include Kahana-like momentum or energy depen- 
dence to describe the results of positron measurements. 
More recently, theoretical prescriptions for the enhance- 



ment have been developed that represent a significant 
departure from the homogeneous electron gas or local- 
density approaches, based on, for example, the general- 
ized gradient approximation [5^, Bloch- modified ladder 
^ or weighted-density approximation [19 . 

Owing to the different screening properties of d and 
s-p electrons, efforts to include a character, or state- 
dependent enhancement function have been applied to 
several transition metals and their alloys [181 I2QH22] . 
Sob applied such a scheme to data measured on a poly- 
crystalline FeAl alloy, finding a de-enhancement of the d 
states by a factor of ~ 2.2 compared with the s-p states 
[18 , whereas the application of the same procedure by 
Svoboda and Sob [20] to CuZn was found to favor a re- 
duction by a factor of ~ 1.5. Theoretically, such explicit 
state-dependence is rarely included, although for fiat d- 
bands it is implicitly present in any energy-dependent 
model. Recently, Barbiellini and co-workers have devel- 
oped a theoretical and ab initio state-dependent prescrip- 
tion for calculating the enhancement in a general system 
[5", '6 , which is based on the state-dependent annihilation 
rates calculated within the generalized gradient approxi- 
mation (GGA). Although it has been demonstrated that 
the effects of enhancement do not shift the location in 
k-space of the Fermi breaks in positron measurements 
[23] , the infiuence of the theoretical treatment of the en- 
hancement, when rigid-band like shifts are applied to the 
electronic structure and compared with experiment, has 
not yet been investigated. 

Here, we tackle the problem of describing the enhance- 
ment of the positron annihilation rates from an experi- 
mental perspective. Employing a state-dependent (SD) 
model for the enhancement similar to that of Ref. [5^ , we 
simultaneously fit both the FS and the enhancement from 
ab initio electronic structure calculations to positron data 
directly in k-space in order to obtain a quantitative mea- 
surement of the enhancement in metals. Additional com- 
parisons with the calculational scheme of Ref. [5^ are used 
to quantitatively assess the applicability of such an SD 
enhancement model for electron-positron momentum dis- 
tributions. In particular, the accuracy of rigid-band-like 
approaches in obtaining more realistic representations of 
the experimental FS are found to be sensitively depen- 
dent on the particular enhancement employed in the cal- 
culation. 

The organisation of this paper is as follows. In Section 
[TTj we introduce the method employed in this paper. In 
Section |III[ we apply this fitting technique to some 3d 
transition and noble metals, namely V, Cr and Ag, and 
in Section |IV| we address the simple metal Al. Finally, 
in Section [V[we apply and investigate a correction to the 
existing theory of Ref. [5 that provides useful predictive 
power as a general model of enhancement in both transi- 
tion metals as well as simple metals. The application of 
this correction to Mo is shown to quantitatively explain 
the difference in the momentum distributions of Cr and 
iso-electronic Mo that is observed despite the similarity 
in their FS. 



II. METHOD 

A. State-dependent enhancement 

The quantity measured by two dimensional angular 
correlation of (electron-positron) annihilation radiation 
(2D-ACAR) experiments is a once-integrated projection 
(along a suitable crystallographic direction) of the so- 
called two-photon momentum density, p^^(p). 



P'^ip) 



Y,n^\C^,G\'S{p-k-G), (1) 



where n^ is the electron occupation density of state 
i = { j, k} {j is the band index), CiG are the coefficients 
of a plane- wave expansion of the product of the electron 
and positron wave functions, in which G is a vector of the 
reciprocal lattice, and the ^-function expresses the con- 
servation of crystal momentum. In a 2D- AC AR measure- 
ment, the 3D quantity expressed in Eq. [l] is integrated 
along a particular direction to yield a 2D projection of 
p^^(p), and the projected axis is usually chosen to be 
a suitable high-symmetry crystallographic axis. The FS 
enters Eq. [l] as discontinuous breaks in the momentum 
density when p traverses i occupied (n^ = N) to i unoc- 
cupied {ui = N — 1) (i.e. when the band crosses the Fermi 
energy). The folding of crystallographically equivalent p- 
points of momentum using the so-called Lock-Crisp- West 
procedure [26 yields the 'reduced momentum density' 
(RMD), p27(k), 



p^-(k) = ^|Q,cp. 



(2) 



i,G 



The Ci^G of Eq.fT]can be written in terms of the single- 
particle electron and positron wave functions, '0i(r) and 
'0+(r) as, 

Ci,G = /d^rexphKk + G) •r]Vr(r,r). (3) 

Here, '^^^(r, r') is the electron-positron pair wavefunction 
for state z. 



C''(r,r')=^i(r')^+(r)vS(r), 



(4) 



where 7^ (r) is the state-dependent positron enhancement 
factor (for the state i). Setting 7 = 1 in Eq. His equiva- 
lent to the independent particle model (IPM), although 
it should be noted that the effects of the positron wave- 
function are still included in that case. The usual pa- 
rameterizations of the enhancement, for example the BN 
or JS models, are based on local-density parameteriza- 
tions, in which 7i(r) = 7(r) is a function only of the 
unscreened local electron density at the location of the 
positron. Other state-dependent prescriptions exist (e.g. 
^ITnlSj), although these have relied on the empirical de- 
termination of the state dependence of the enhancement. 
Barbiellini and co-workers [5 have proposed a theoret- 
ical prescription for applying a state-dependent positron 



enhancement factor to ab initio calculations of the elec- 
tronic structure and momentum density. In their scheme, 
ji is obtained through the partial annihilation rates, such 
that. 



7i = Xi/X. 



IPM 



(5) 



where A^ is the partial annihilation rate of the state i 
including correlation effects, and X]^^ is the partial an- 
nihilation rate due to the IPM. The total annihilation 
rate. A, may be calculated from (here shown for the local 
density approximation, LDA), 



A 



I / d^r n+(r)n(r)7(r). 



(6) 



where r^ is the classical electron radius, c is the speed 
of light and n+(r) is the positron density. In their cal- 
culations, the GGA was used for the calculation of A^, 
which successfully reproduces the experimental annihila- 
tion rates rather well |TOl [27| . 



In this way, 7/ is a universal quantity, representing the 
partial enhancement of a state with character /. The de- 
gree to which it is enhanced depends on the coefficients 
of the wavefunctions in the LMTO calculation. The en- 
hancement, then, of a pure state of atomic species n and 
orbital character / is given by 7^,^ = 7n •7^ Note that the 
band characters (atomic species and orbital character) 
are strongly k-dependent, and of course vary from band 
to band due to hybridization with other states, so our en- 
hancement model is a general state-dependent model for 
the enhancement (see, for example. Fig. fl]), but has its 
origins in the convenient properties of the LMTO wave- 
functions. 

The contribution due to core annihilations is an impor- 
tant consideration for positron lifetime calculations [28]. 
However, the core contribution is small and relatively in- 
dependent of k across the first Brillouin zone (BZ), and 
can safely be omitted from this calculation. Instead, the 
contribution from core states in the data is described by a 
uniform background in the subsequent fitting procedure. 



B. Practical approach 



C. Minimization procedure 



We begin by computing the ab initio electronic struc- 
ture using the linearized muffin-tin orbital (LMTO) 
method, within the atomic sphere approximation and 
including combined correction terms [9j. The Ci^o co- 
efficients of Eq. [3] are then computed within the IPM 
(equivalent to setting 7 = 1 in Eq. [4|, unfolded in such a 
way as to resolve the individual contribution owing to the 
atom index (n), and orbital angular momentum quantum 
number (/), 



^IPM 



E^IPM 
n,l 



(7) 



The momentum density in the first Brillouin zone (i.e. 
the RMD, Eq. |2| is computed for the IPM by, 



pf^ik) = constant x ^ \{J2cl^n,i)\' 



(8) 



n,l 



For the enhancement, we introduce the quantities jn,i 
that describe the enhancement of a state of atomic 
species n and of orbital angular momentum / (/ = 
s^p^d^f). These can then be incorporated into the cal- 
culation of the RMD by. 



p^^(k) = constant x ^ \(J^ 

G n,l 



(9) 



Note that in the above equation, the ji multiply the 
Ci^G,n,i coefficients, which are inside the sum over in- 
cluded G- vectors, and so the RMD must be re-computed 
for each 7/ and cannot be expanded into a sum of con- 
tributions to the momentum density from different l- 
orbitals. 



In the rigid-band approach, the agreement between ex- 
periment and theory is iteratively maximized with re- 
spect to a rigid shift of one or more of the energy bands 
(typically those that constitute the FS), until conver- 
gence at the minimum of the goodness-of-fit parameter is 
achieved. This is similar to the method of Ref. [T4], how- 
ever, there are some important differences. In Ref. [ll], 
the radial anisotropy of the two-photon momentum den- 
sity in p-space served as the comparative quantity, and 
the enhancement was fixed to that chosen in its ini- 
tial calculation (in that case, the JS model was used). 
Here, we perform our comparison in k-space, correspond- 
ing to the Lock- Crisp- West-folded data, and explicitly 
include enhancement of the form outlined above (SD 
model) in the fitting. The advantage of operating in 
k-space (aside from the smaller array sizes involved) is 
principally that we are sensitive directly to the projected 
Fermi breaks, rather than the many weaker FS signa- 
tures that are distributed throughout the p-space spec- 
trum. An additional consideration, however, is the con- 
tribution from higher momentum components (Umklapp 
processes), whose enhancement has presented a challenge 
for theoretical models (see, for example, Ref. [7 ). It is 
noted that operating in k-space involves the folding of the 
Umklapp contributions into the first BZ, both experimen- 
tally and theoretically, and that any non-trivial behav- 
ior of these contributions is subsequently lost. However, 
we have checked our results with the equivalent p-space 
spectra, and in particular near the Umklapp regions (as 
well as its integral, which represents an analogue of the 
coincidence Doppler broadening spectra). Crucially, we 
find the data are equally well-reproduced using such a k- 
space approach as they are with the traditional JS model. 

The fitting parameters constitute the energy shift, 5j^ 



for each band in the fit (typically those that cross E'p), 
two scaling parameters for each experimental projec- 
tion, A^ and Sm, that approximately relate to the k- 
independent core contribution and the number of counts 
in the 2D- AGAR spectra respectively, and the enhance- 
ment parameters, 7^,^, of which there are typically three 
for simple systems. These are simultaneously adjusted 
using the MINUIT package [29 and the computed k-space 
density is compared with the data until convergence is 
reached. Note that we fit the ratios of the enhancement 
parameters, absorbing their magnitude into the scaling 
parameters (the absolute magnitude of the enhancement 
parameters is indistinguishable from the scaling parame- 
ters in the data). 

Our definition and treatment of the scaling parameters 
have an important consequence. As mentioned above, we 
do not treat the enhancement of the core electrons, pre- 
ferring to concentrate on the description of the shape of 
the RMD. Such an approach means that good agreement 
can be obtained with the IPM if we consider a nega- 
tive contribution from core annihilations. Whilst this is 
clearly unphysical, it stems from the strong overestima- 
tion of the enhancement of deeply-bound electrons within 
the IPM. Here, we are most interested in the band prop- 
erties of the momentum distribution (i.e. its shape) ^ and 
in particular its FS signatures. We point out that in 
the following discussion, even when the IPM appears to 
give reasonable agreement with our data, the agreement 
with positron lifetime measurements (see, for example, 
Refs. [13j 28j 30 ) would be very poor, in contrast to the 
other enhancement models that are addressed here. 



III. TRANSITION AND NOBLE METALS 

The transition metals and their alloys have tradition- 
ally been the subject of the bulk of experimental investi- 
gations of the FS, and a good description of the electron- 
positron momentum density and enhancement of such 
systems has been vital in understanding 2D- AGAR, and 
indeed coincidence Doppler broadening [6] [31] [32] data. 
The JS model was specifically developed with transition 
metals in mind, and is generally thought to provide a 
good description of the enhancement for d-electron den- 
sities (with rs ^ 1.8) [T2HII]. Here, we begin by apply- 
ing the SD enhancement model described above to some 
metals whose FSs have been accurately determined via 
quantum-oscillatory methods (Ag, V) and one whose FS 
is inaccessible to conventional FS probes (paramagnetic 
Gr), first concentrating on the "raw" LMTO calculations 
of the RMD. Gomparisons are made with both the IPM 
and the JS model for enhancement, as well as a simpli- 
fied version of the Barbiellini-Alatalo [5, 6 enhancement 
scheme. Following this, we rigidly fit the bands to the 
experimental data to obtain an experimental measure- 
ment of the FS, in order to assess the sensitivity of this 
approach to the FS. 

Three 2D- AGAR projections ([100], [110] and [111]) 
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FIG. 1: (color online) The state-dependent enhancement of 
Ag from our model, shown for two energy bands along the 
path F-X in the BZ. (a) The dispersion of bands 1 and 6; 

(b) the enhancement from our fit to the experimental data; 

(c) and (d) the character of bands 6 and 1 respectively. Note 
that band 6 crosses ^f, shown by the dotted line in (a); above 
which the enhancement is unphysical. 



were obtained from a single crystal of Ag at ~ 70 K, with 
a resolution full width at half maximum of 0.71 x 1.11 
mrad in the Px and Py data axes respectively (correspond- 
ing to - 12% X 19 % of the BZ of Ag). For V, four pro- 
jections were obtained along the [100], [110], [210] and 
[211] directions at room temperature with a resolution 
of 1.11 X 1.33 mrad (with the exception of the [110] di- 
rection, which was collected at ~ 24 K with resolution 
0.83 X 1.11 mrad). Paramagnetic Gr was measured along 
the [100] and [110] directions at 353 K, well above the 
Neel temperature (with a resolution function the same as 
the room temperature V measurements). For each sam- 
ple presented in this manuscript, the 2D- AGAR spectra 
have been carefully checked to confirm the absence of any 
defect or impurity signatures in the spectra. 

LMTO calculations were performed over 1505 k-points 
in the irreducible wedge of the face-centered cubic BZ 
for Ag, and over 6201 k-points in the irreducible wedge 
of the body-centered cubic BZ for V and Gr. For each 
material, the RMD was computed for both the IPM and 
the JS parameterizations of the enhancement, and convo- 
luted with the appropriate experimental resolution func- 
tion. This was compared to the experimental data with 
adjustments only to the scaling parameters, Sm and A^. 
The SD model of the enhancement was then obtained 
by simultaneously fitting the orbital enhancement fac- 
tors, ji. In each case, a goodness-of-fit parameter, x?ed' 
was computed as a weighted average of that from each 
experimental and theoretical projection. Gomparisons 
were also made for the raw band calculation with the 
model of Refs. [5, 6 by computing the annihilation rates 
associated with each orbital both within the GGA [lOj 
and IPM, where the enhancement of each orbital in this 
model is the ratio of these annihilation rates. However, it 
should be noted that this scheme still invokes a parame- 
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FIG. 2: (color online) Comparison between experimental k-space momentum density for (a) Ag [100] projection, (b) V [110] 
projection and (c) Cr [110] projection and the computed raw band calculation of the RMD for IPM, JS and SD models. High 
symmetry points in projection have been labelled. Note that in (a) the horizontal and vertical axes are (110) crystallographic 
axes; the (100) axes are along the diagonal. 



terization of the enhancement in terms of the electron-gas 
parameter. To compare with our experimental values, 
we integrate over all k-states to obtain the annihilation 
rate from all electrons of orbital quantum number /. In 
Eq. [6J this corresponds to substituting the (partial) elec- 
tron density due to each 5, p, d and / orbital, n/(r) for 
n(r) to obtain Xi in the GGA, rather than the method 
of Ref. [5J in which nj,k(r) is used. In the following, we 
refer to this as the simplified Barbiellini-Alatalo (SBA) 
model (in which the 'simplified' reflects the integration 
over all k-states). 



A. Raw band calculations 

For the raw band calculations, the comparison between 
theory and experiment depend on (i) a good ah initio de- 
scription of the electronic structure, and in particular 
the FS, and (ii) a reliable understanding of the electron- 
positron enhancement factor. The well-known FS of Ag 
consists of just a single sheet that is only slightly per- 
turbed from the free-electron sphere, most notably along 
the [111] direction where the FS intersects with the BZ 
boundary to form a neck at the L-point of the BZ. Band 
structure calculations within the LDA reproduce the pre- 
cise measurements of quantum oscillations [33H35^ very 
well and it therefore provides an excellent candidate in 
which to test models for the enhancement, for which 
the JS model would be an obvious choice. As demon- 
strated by Fig. [2k and [Sk, the experimental data are 



well-described by the raw band calculations of the RMD, 
in which even the IPM (including only positron wave- 
function effects) works reasonably well. Quantitatively, 
as demonstrated in Table [ij the JS enhancement model 
is found to improve the agreement between experiment 
and theory, particularly near the projected TX and L 
points of the BZ. However, the current SD model is able 
to bring the theoretical RMD into much closer agreement 
with the data by de-enhancing the s and d states relative 
to the p states. For Ag, the bands below E^ are pre- 
dominantly (i-bands, with some hybridization with the 
5s state, but the band that crosses Ep has substantial 
p character. The de-enhancement of d states relative to 
the sp bands is well-known [12] and is attributed to the 
relative localization of d electrons, particularly near the 
top of the d-bands. That the p enhancement appears 
to be quite strong can be explained by a Kahana-like 
momentum enhancement, in which those electron states 
nearest kp are most enhanced. In Fig. [l] the measured 
SD enhancement is plotted along F-X, accompanied by 
the band dispersion and character. As can be seen from 
Fig. [TJd, the enhancement of band 6 grows substantially 
as the band approaches the Fermi level, replicating the 
Kahana-like momentum dependence of the enhancement. 
This is captured in our model by the enhancement of the 
p-like states; as demonstrated by Fig. [T]^, the enhance- 
ment of both bands 1 and 6 closely follow their respective 
p character. Note that, owing to the weak contribution 
from high- lying / states, a good quantification of their 
enhancement is not possible. 
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TABLE I: The results of the fit between the different pa- 
rameterizations of the enhancement and the data. For the 
SD model, the ji for each state obtained from the fit is also 
given, normalized to 7p = 1. The 7^ for the SBA model are 
determined from the partial annihilation rates described in 
Eq. [5] The errors in the fit of the ji of the SD model are 
^ zbO.Ol. Note that the higher statistical precision of the V 
data yields a relatively large x?ed parameter. 



The situation is more complicated for V and Cr, for 
which the details of the near-£^F electronic structure, in- 
cluding the precise dimensions of the FS, are either not 
well-reproduced by our band calculations (V), or have not 
been accurately determined (paramagnetic Cr). Figs.|2b 
and [2J3 show representative 2D-ACAR projections of V 
and Cr (along the [110] direction) respectively, compared 
with the corresponding theoretical quantities. The FS of 

V is composed of two sheets. The first sheet (originat- 
ing from band 2 of predominantly 3d character) forms 
a small F-centred hole octahedron that encloses ~ 0.12 
holes and has remained unobserved in quantum oscilla- 
tion data, although its presence has been confirmed by 
2D-ACAR measurements [36 - 38^ (this sheet is visible in 
the experimental data in Fig. ^ at the projected F7V 
point). Band 3, on the other hand, experiences apprecia- 
ble hybridization with the 4p states above E^ and forms 
a F-centred jungle-gym hole FS as well as some hole el- 
lipsoids that are centred at N. These TV-hole ellipsoids 
can be clearly seen in the data of Fig. [2)3 at the projected 
A^-point of the BZ, where the density experiences a local 
dip due to their presence. These features in the IPM and 
JS calculations are predicted to be substantially too large 
and too strong (see Fig.lsb), and the enhancement of the 
high-density surrounding region is not well reproduced. 
Although, as pointed out by Jarlborg and Singh, the en- 
hancement is expected to be less important for a less-full 
(i-band, the JS enhancement is actually found to perform 
worse than the IPM for V (see Table ^, at least in the 
shape of the distribution (positron lifetime predictions 
are substantially better described by the JS model [13 ). 

V was used as a test material by Jarlborg and Singh in 
their presentation of the JS enhancement model, in which 
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FIG. 3: (color online) The raw-band RMD of (a) Ag, (b) V 
and (c) Cr shown in Fig. [2] shown here along a path in the 
BZ. 



they comment that the IPM already provides a reason- 
able description of the momentum density [12] , and that 
their model offered only weak improvement. However, 
their comparisons were made with electronic structure 
calculations where the bands had been rigidly shifted to 
agree with de Haas- van Alphen (dHvA) measurements of 



the A^-hole ellipsoids [39^. Indeed, in Fig.jSJD it is obvious 
that the dimensions of the dips in the momentum density 
along NH-N-HN are incorrectly placed with respect to 
the data. In Section IIIB[ we address such inconsisten- 
cies by rigidly shifting the bands to improve agreement 
between experiment and theory. 

Cr neighbors V in the periodic table, having an ex- 
tra electron, yet its paramagnetic FS has remained rel- 
atively unexplored experimentally, principally owing to 
the emergence of an ordered spin-density wave phase be- 
low ~ 312 K [40 , where the high temperature precludes 
quantum oscillatory measurements in the paramagnetic 
phase, and strong spin fluctuations appear to suppress 
the measurement of the nested sheets in the ordered 
phase [41]. Theoretically, the FS is composed of three 
sheets, the first of which (band 3) contributes some small 
electron 'lenses' midway between F and H. Band 4 forms 
some N hole ellipsoids and ilf-centred octahedra, whereas 
in band 5 there is a F-centred electron 'jack'. Molybde- 
num, isoelectronic to Cr, shares a similar FS topology, 
in which the A/"- hole ellipsoids and the electron jack can 
be clearly visualized in the [110] -projected k-space den- 
sity of 2D- AGAR measurements |42l [43]. For Cr, these 
features, shown in Fig. [2J3 near the projected N points, 
are obscured in the measurement, presumably owing to 
enhancement effects [44]. Indeed, the A^-hole ellipsoids 
are more evident in the IPM and JS projected densities 
than they are in the data. Overall, the agreement be- 
tween the IPM and JS calculations of the RMD and the 
data is reasonable (see Table |l| , but is particularly poor 
near the projected A^-points of the BZ as well as midway 
along the NH-TN path (see Fig. [sj^). Here, the knobs of 
the electron jack project on top of one another, and the 
IPM and JS do not predict the de-enhancement of the 
momentum density very well in this part of the BZ. 

The application of the SD model, however, consider- 
ably improves the agreement between experiment and 
theory by substantially de-enhancing the s and d states. 
In Figs. [3J3 and [3J3, this can be most clearly seen at the 
projected A^-points of the BZ, as well as the momentum 
density near the TN points. For both V and Cr, the hy- 
bridization of the valence states with the unoccupied 4p 
states indicates the importance of the p electrons in de- 
ciding the topology of the FS, and their proximity to E^ 
means they have a strong impact on the enhancement of 
the momentum density in 2D-ACAR measurements. Pre- 
vious non- iterative comparisons of orbital- weighted band 
theory and 2D-ACAR data have been made in p-space 
for V [22], in which a de-enhancement of the s and d 
states by ^ 0.8 relative to the p states was favored. Our 
results are close to these, where we obtain 0.69 and 0.78 
for s and d states respectively, corresponding to a slightly 
greater de-enhancement of the lower-lying s states. Sim- 
ilarly for paramagnetic Cr, Matsumoto and Wakoh [21] 
estimated (also non-iteratively) that the Cr d states were 
de-enhanced by ~ 0.67 relative to the sp states (which 
were considered together). Our results correspond well 
with their findings, where we obtain 0.82 and 0.54 for s 



and d states, relative to the p states. The stronger en- 
hancement of the Cr s states (compared with V) may be 
explained by the higher occupation of the 4s states in Cr 
(they are almost twice as occupied in Cr). 

Finally, we comment on the predictions for state en- 
hancement made by the SBA model. Apart from Ag, 
which has a much higher d electron density than either V 
or Cr, the orbital enhancement ratios are predicted to be 
very similar (see Table [l|), with a weak de-enhancement of 
the 5-states and a modest de-enhancement of the (i-states 
relative to the p-states. For Ag, a rather more exagger- 
ated de-enhancement is predicted for the (i-states. Qual- 
itatively, these results are in agreement with our mea- 
sured values, but differ substantially in magnitude and 
lead to a slightly higher Xred parameter than the current 
SD model. Nevertheless, the SBA model provides better 
agreement with the data than either the IPM or the JS 
model for Ag and Cr, supplying a more robust predictive 
scheme for computing the enhancement in 2D-ACAR mo- 
mentum distributions. That it does not fair so well for V 
is mostly accounted for by the rather larger corrections to 
the LDA band structure that are required for V, a topic 
that will be returned to in the next section. Here, we em- 
phasize that the SBA model is expected to improve the 
agreement with the data over the IPM or JS models when 
extensions to the LDA, such as non-local potentials [45] 
or self-energy corrections [46^, that improve the descrip- 
tion of the FS are included. The most probable origin for 
the discrepancy between the SBA scheme and our mea- 
sured SD enhancement is the omission of a Kahana-like 
energy-dependence or momentum-dependence in the pre- 
dictive scheme of Ref. [5 . As already highlighted, some of 
the results for the enhancement in the current SD model, 
particularly the apparent strong p enhancement in V and 
Cr, reflect a Kahana-like enhancement of those electron 
states near £^f- In our model, where a state is not too dis- 
persive in energy, this is naturally captured by enhance- 
ment of that state. The authors of Ref. [5] comment 
that the Kahana-like energy- or momentum-dependence 
appears to be less important than the state-dependence 
from their results, a conclusion that this work substanti- 
ates, but these results suggest that including such en- 
hancement could produce a good improvement in the 
agreement between experiment and theory. In Section 
\V\ we apply just such an energy-dependent term to the 
SBA enhancement factor, and demonstrate the improved 
predictive capacity of such a model. 



B. Rigid-band fit 

In 2D-ACAR investigations of the FS, the traditional 
method of extracting the FS from experimental data is 
to contour the data at a level that corresponds to ex- 
trema in the first derivative of the data, and it is well 
known that enhancement effects do not shift the loca- 
tion of these breaks [23]. While first-principles calcula- 
tions are often able to make excellent qualitative predic- 




wx 




FIG. 4: (color online) Comparison between experimental data 
for (a) V [110] projection and (b) Cr [110] projection and the 
rigid band fit to the RMD for the SD model, shown in the 
same way as Fig. [2] 



tions about the nature of the Fermi surface, when sub- 
ject to detailed scrutiny in light of precise experimental 
data it is often found that quantitative differences exist. 
Shortcomings in the approximations used in the calcu- 
lations (e.g. exchange-correlation functional, neglect of 
relativistic effects) mean that in reality it is difficult to 
get the Fermi surface correct. These differences can of- 
ten be reduced or eliminated by small shifts of the rel- 
evant bands with respect to the Fermi level, and it has 
recently become feasible, and indeed quite common, to 
'tune' a band-theoretical calculation in this way (e.g. see 
Refs. [T4||24l [25] ) . Such an approach requires an accurate 
description of the positron enhancement, if conclusions 
regarding details of the FS itself are to be drawn from 
such a fit, and we now turn our attention to investigate 
the behavior of our SD model applied to such detailed 
FS studies. 

For the rigid-band fit of the electronic structure, the 
Fermi level for each band near Ep was fitted to the data. 
In the case of the SD enhancement model, the orbital en- 
hancement factors were fitted simultaneously. The x?ed 
was computed as before, and the number of electrons en- 
closed by the fitted FS (i.e. the occupied fraction of the 
Brillouin zone) was obtained. 

The results of the rigid-band fit to the data are dis- 
played in Fig. [4J and demonstrate substantial improve- 
ment over the corresponding raw band calculations of 
Fig. p] In Table [ll{ the x?ed ^^ shown for each fit, along 
with the fitted orbital enhancement factors for the SD 
model. Beginning with some general comments, we note 
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FIG. 5: (color onUne) The rigid-band fit to the RMD of (a) 
Ag, (b) V and (c) Cr shown in Fig. HI shown here along a 
path in the BZ. 



that the orbital enhancement parameters obtained in the 
SD model are only moderately adapted as a consequence 
of including the bands in the fit, and the same general 
trends are observed. Additionally, it is noteworthy that 
in al most every case the shift in the energy band (see 
Table III ) is found to be smallest for the SD model (with 
the exception of band 3 of V). 





7s Ip Id 7/ 


X?ed 


Ag 


IPM 
JS 
SD 


0.88 1.00 0.85 (0.76) 


5.53 
4.84 
4.09 


V 


IPM 
JS 
SD 


0.61 1.00 0.63 (0.49) 


7.77 
8.40 
3.79 


Cr 


IPM 
JS 
SD 


0.83 1.00 0.61 (1.27) 


2.52 
1.85 
1.28 



TABLE II: The results of the rigid-band fit between the dif- 
ferent parameterizations of the enhancement and the data, 
presented in the same way as Table [l] The shifts in the en- 
ergy bands for each fit are shown in Table [ill] 



Fig. [5] shows the RMD along the same path through 
the BZ as in Fig. [3J and we will now concentrate in more 
detail on the agreement between experiment and theory. 
As indicated by the small change in electron count of the 
shifted bands, the change in the FS itself is small, owing 
to the appreciable dispersion of band 6 at £^f, and the 
shift in the Fermi wavevector is A/cp ~ 0.02 (27r/a) (just 
~ 15% of the resolution function). 

The improvement is much more dramatic for V (Fig. 
[5b), in which the size of the A/"- hole features in the data 
is now well-described by all of the enhancement models, 
stemming from opposite shifts in bands 2 and 3. After 
rigidly shifting the bands, the IPM and JS demonstrate 
similar shifts of the energy bands and a similar goodness- 
of-fit parameter, leading to an excess in occupied volume 
of 0.14 and 0.10 of an electron The improvement in the 
SD model is more pronounced, however. The d bands 
are well known to be placed too low by the LDA with 
respect to sp bands. As noted in Ref. [14], band 2 (of 
predominantly d character) is pushed up in energy by 
the fit towards the higher p character of band 3, which is 
pulled down by the fit, correcting this tendency. 

The energy shifts of the bands are also in good agree- 
ment with quantum oscillations. Comparing the semi- 
axis radii of the ellipsoids with high-quality dHvA pa- 
rameterizations (see Table IV), we find much improved 
correspondence with experiment than the raw calcula- 
tion, and indeed they compare favorably with the shifts 
of Ref. ^14^. The differences in the results of Ref. [14 and 
the current k-space approach reflect the different sensi- 
tivity of the two techniques to specific features of the data 
(for example, compare the N-P radius with the N-H ra- 
dius). It is also worth mentioning that the jungle-gym FS 
also originates from band 3, and that this will also con- 
tribute to the shifts of this band, and so considering the 
orbits about the A/'-ellipsoids alone may be misleading. 
Unfortunately, there is a dearth of data for this sheet of 
FS, and comparisons are hard to draw. As a final point, 
quantum oscillations appear to be relatively insensitive 
to the FS of band 2, whereas we find a strong dependence 





band shifts / mRy 


electron +/- 


Ag 




band 6 




IPM 
JS 
SD 


-21.2 
-18.3 
-14.2 


+0.08 
+0.07 
+0.05 


V 




band 2 band 3 




IPM 
JS 
SD 


+22.3 -15.8 
+26.5 -15.4 
+ 18.8 -16.4 


+0.14 
+0.10 
+0.17 


Cr 




band 3 band 4 band 5 




IPM 
JS 
SD 


-25.1 +15.6 +13.2 
-22.9 +12.2 +13.2 
-18.5 +5.0 +13.2 


-0.02 
-0.01 
0.02 



TABLE III: The shifts in the energy bands for each of the 
rigid-band fits. Also shown is the change in electron count 
in the BZ due to the fit. Note that for Cr band 5 (that just 
grazes ^f) is completely expelled by all of the fits. The errors 
in the shifts of the bands are in each case < 1 mRy. 



of our fit to that band, in agreement with other positron 
studies in V [36, 37 . 

Finally, for Cr (Fig. Isfc), the data are reasonably well 
described by the SD model throughout the BZ, whereas 
the IPM and JS models struggle near the FA^ points in 
both the raw band calculations and the rigid-band fits. In 
the absence of high-precision FS data for paramagnetic 
Cr, owing to the ordering temperature (Tn ~ 312K) of 
the spin-density wave, a robust comparison can instead 
be made of the nesting vector of the paramagnetic FS 
that is widely believed to determine the ordering (and 
has remained difficult to establish experimentally). High- 
resolution neutron diffraction measurements have estab- 
lished the ordering vector to be Q = (0, 0, 0.9516) (27r/a) 
(see Ref. [40^ and references therein). Our raw LMTO 
calculations predict (via a computation of the static sus- 
ceptibility, Xo(q), see for example Ref. [47]) a nesting vec- 
tor q ~ 0.930 (27r/a) (see Fig. [6]), which is rather smaller 
than the neutron measurements. Since quantum oscilla- 
tions are precluded in the paramagnetic phase (and have 
recently remained unobserved from the relevant FS sheets 
in the ordered phase owing to strong spin-fluctuation in- 
duced scattering [41 ), the only data on the FS that has 
been capable of extracting this nesting vector have been 
some recent angle-resolved photoemission measurements 
on Cr(llO) thin films [48", "49^, in which a nesting vector 
of q ~ 0.950 ± 0.005 (27r/a) is reported, in very good 
agreement with neutron measurements. Our rigid- fit to 
the data (Fig. [6| culminates in a FS nesting vector of 
q ~ 0.950 ± 0.002 (27r/a), where the error quoted is the 
combined error from shifting the two bands to match the 
experimental results, representing the highest-precision 
experimental confirmation of the relevant dimensions of 
the FS of paramagnetic Cr from a bulk measurement. 
In contrast to this excellent agreement, the shifts of 
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direction 


LMTO 


Ref. [39] dHvA 


SDfit 


Ref. [14] fit 


N-P 

N-r 

N-H 


0.257 
0.254 
0.168 


0.223 
0.212 
0.176 


0.224 ±0.002 
0.204 ±0.002 
0.146 ±0.001 


0.245 
0.231 
0.160 



0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96 

q,oo(27r/a) 



FIG. 6: (color online) The real (top) and imaginary (bottom) 
parts of the static susceptibility, Xo(q), of paramagnetic Cr, 
calculated for the raw band LMTO calculations as well as 
the results of the rigid-band fit to the data with the IPM, 
JS and SD models of enhancement. The dashed vertical line 
represents the peak in the real part of the susceptibility. The 
inset shows a slice of the FS through the (001) plane, with 
the arrow depicting the nesting that gives rise to the peak in 
3m Xo(q) between the hole (outer, red) and electron (inner, 
blue) FS sheets. 



the bands obtained using the IPM and JS models sug- 
gest nesting vectors of q ~ 0.924 ± 0.002 (27r/a) and 
q - 0.932 ± 0.002 (27r/a) respectively. 

The conclusions we draw from this section are the fol- 
lowing. First, our approach provides a robust empiri- 
cal means of measuring the orbital electron-positron en- 
hancement factors, that are truly state-dependent (i.e. k- 
dependent). Second, this measurement is not so strongly 
dependent on the accuracy of the band calculation, being 
rather more sensitive to the overall shape of the momen- 
tum distribution. Third, simultaneous fitting of the en- 
ergy bands and the orbital enhancement lead to a tuned 
FS that is in better agreement with other FS data than 
is the raw band calculation, as well as in good agreement 
with previous p-space fitting approaches. Finally, the 
band shifts that are required to reproduce experimental 
data (that is in better agreement) are generally smaller 
for the SD model than the other approaches investigated 
here, indicating that artificially large rigid shifts in the 



TABLE IV: Comparison of the semi-axis radii (in units of 
271 /a) of the N hole ellipsoids from our raw LMTO calculation 
and the fitted SD momentum density for V. Comparisons are 
made with the high-precision parameterizations of dHvA data 
of Ref. [39] as well with the fitting technique (also applied to 
2D-ACAR data) employed by Ref. Il4j. The errors reflect 
the error in locating the minimum of the fit with respect to 
the shift in the bands. Note that the dHvA radii rely on the 
assumption of perfect ellipsoids. 



rigid fit 




FIG. 7: (color online) Comparison between experimental data 
for Al projected down the [111] axis and the computed RMD 
for IPM, JS and SD models, (a) The raw LMTO band cal- 
culation, and (b) the results of the rigid fitting of the energy 
bands. The high symmetry points in projection (FL, X and 
W) are shown in (a), and the boundary of the first BZ is 
marked by the dotted line. 



bands can develop as a consequence of an inadequate de- 
scription of the enhancement. 



IV. SIMPLE METALS 

We now turn to the other regime of enhancement, in 
which the bands, of sp character, are closer to the nearly 
free electron model. Aluminium and the alkali metals 
(and their alloys) provide a more stringent test for the 
SD model. The electron-gas parameter of Al is r^ = 2.65, 
above the point at which the JS (which does not conserve 
the low-density limit) and the BN (which does) begin to 
diverge; as a consequence JS is not expected to perform 
well here. 



11 





7^ 7p Id 7/ 


■A/2 
Xred 


raw band 


IPM 


_ 


15.40 


JS 


- 


17.20 


SD 


0.53 1.00 (1.29) (0.60) 


2.29 


rigid fit 


IPM 


_ 


3.93 


JS 


_ 


3.92 


SD 


0.60 1.00 (1.09) (0.72) 


2.12 



TABLE V: The results of the fit between the different param- 
eterizations of the enhancement and the data for Al, shown 
in the same way as Table [III The band-shifts that accompany 
the rigid-band fit are shown in Table IVl] 



The FS of Al is composed of two bands, one which 
forms a F-centred hole sheet from band 2 that lies com- 
pletely in the first BZ, and from band 3 a so-called 'dis- 
membered monster', that consists of square electron rings 
that run the length of the edges of the first BZ except 
at the corners (VF-points) [50]. Two 2D- AGAR projec- 
tions along the [110] and [111] directions were measured 
at room temperature and compared with LMTO calcula- 
tions performed over 1505 k-points in the irreducible BZ 
using the IPM, JS and SD enhancement as before. 

In Fig. [tI the data for the [111] projection is shown 
alongside the LMTO calculations of the RMD. The FS 
structure can be clearly seen in the data (shown in the 
bottom right panel), where the low density in the cen- 
ter reflects the band 2 hole sheet, and the higher density 
at the edges of the projected BZ come from the electron 
rings of band 3. At the corner of the projected BZ (near 
the FF-point), the particularly high region is due to the 
projection of the rings in neighbouring zones along the 
(111) directions. As can be seen in the left panels of Fig. 
[T) the IPM and JS models are particularly poor at de- 
scribing the enhancement at the edges of the zone that 
connect these strong features. Moreover, a small local 
peak at FL that is predicted by both IPM and JS is not 
observed at all in the data. Quantitatively, as might be 
expected from the electron density of Al, the JS model 
fairs poorly for the raw band calculation, and even worse 
than the IPM (see Table [V]). The SD model, however, 
does an excellent job of describing the RMD of Al, cor- 
rectly accounting for the absence of the local peak at FL 
and the connectivity of the strong features near FF, and 
leading to an almost order-of-magnitude improvement in 
the x?ed parameter. Here, the s-states are de-enhanced 
substantially, presumably owing to them lying very low 
in energy. Unlike the previous d-electron systems, the de- 
enhancement of the (unoccupied) d states is not observed 
for Al. 

When the bands are fitted, the agreement between 
data and theory for each model is very good. However, 
for the IPM and JS models the local peak at TL persists, 
albeit at a much weaker amplitude. Moreover, consistent 





band shifts / mRy 


electron +/- 




band 2 band 3 




IPM 
JS 
SD 


-39.2 -30.0 
-40.0 -33.7 
-11.9 -10.2 


+0.21 
+0.22 
+0.07 



TABLE VI: The shifts in the energy bands for each of the 
rigid-band fits of Al. 



with the previous conclusions, the shifts in the energy 
bands are substantially larger fo r the IPM and JS mod- 
els than the SD model (see Table VI), leading to an elec- 
tron excess of ~ 0.21 (over a single FS sheet). For the 
SD model, this discrepancy is much reduced, at just 0.07 
electrons. Similarly to Ag, the FS of Al is already well- 
described by the LMTO calculation, and comparisons 
with quantum oscillatory data [51, 52 agree with the raw 
band and SD rigid-band fit to within A/cp ^ 0.03 (27r/a) 
(~ 15% of the resolution function). 



V. PHENOMENOLOGICAL MODEL 

Given the above results, we aim to find a phenomeno- 
logical model that imparts predictive capability on the 
calculation of the RMD. Taking the SD fitted FS as a 
baseline, we attempt to improve on the SB A model of 
the enhancement. The predictions of the SB A enhance- 
ment are, in general, satisfactory, offering a similar de- 
scription of the experimental RMD (in some cases slightly 
better, in others slightly worse) to the JS enhancement 
model. The predictions of the SBA model can be under- 
stood largely from the perspective of the localization of 
the states, in which s and p states experience similar en- 
hancement over the IPM, with the s states in transition 
metals slightly less than p due to their slightly more local- 
ized nature in these systems. The d states are enhanced 
much less in the transition metals, associated with the 
greater localization, and the increasing localization as the 
d-band becomes more filled is reflected by the greater de- 
enhancement of the d-states in Ag when compared with 
either Cr or V. 

The Kahana model for enhancement, applied to a ho- 
mogeneous electron gas and parameterized in terms of 
(A://cf)^, is not expected to work well for (i-band systems, 
in which the effects of the crystal lattice can completely 
hide the Kahana nature. For this reason, Mijnarends 
and Singru (MS) [15] proposed a scheme parameterized 
hy e = {E- £^bot)/(^F - £^bot), where E\,ot is the energy 
at the bottom of the conduction band. 



7 = a + 6e + ce 



(10) 



where a, h and c are constants determined by the electron 
gas parameter Tq . For a parabolic s band this is identical 
to Kahana's formalism. MS demonstrated the applica- 
bility of their prescription for the case of Cu, in which 
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SBA 


SBA-MS 




Xred 


N(Ef) 


b/a 


A/2 
Xred 


V 
Cr 

Ag 


20.41 
2.50 
4.71 


23.80 
9.52 
3.60 


0.700 
0.202 
0.042 


4.30 
1.39 
4.50 


Al 


7.90 


5.04 


0.589 


2.28 



TABLE VII: The linear component of MS-type energy- 
dependent enhancement (b/a) obtained by fitting the SBA 
model to the data. N{Ef) is given in units of states / Ry / 
unit cell, and the quadratic term, 



10 



is set to 0. 



c/a, in Eq. 
The x?ed parameter is given before (SBA) and after (SBA- 
MS) the application of the MS-type enhancement. 



substantial improvement was found (in p-space) with this 
description. The SBA model accounts for the variations 
in enhancement due to the localization of a particular or- 
bital, and its overlap with the positron wavefunction, but 
does not consider the proximity of a state to the Fermi 
level, leading to an over-estimation of the enhancement 
of more tightly-bound, filled s electron shells. Adding 
such a scheme to the SBA model was not found to uni- 
versally explain the variations in enhancement for our 
experimental data without different choices of the con- 
stants a, b and c (in fact, following MS, we choose to set 
a = 1 in Eq. 10 so that b -^ b/a and c 



c/a). 



Of particular interest in this analysis is the enhance- 
ment for V and Cr, which are neighbors in the periodic 
table and would therefore, from the perspective of a ho- 
mogeneous electron gas, be expected to follow similar 
trends in their enhancement owing to their similar elec- 
tron density. As can be seen in Table |TI| the measured 
enhancement of V and Cr is quite different, and yet V and 
Cr can each be well- approximated by a calculation of the 
other's electronic band structure, with a simple extrapo- 
lation of Ef to account for the different band-fillings (i.e. 
the rigid-band approximation works well). The usual pre- 
scriptions for enhancement, in terms of the electron den- 
sity, or even a MS type energy-dependent enhancement, 
fail to predict such different shell enhancements. Sub- 
stituting the measured SD enhancement parameters for 
Cr into the V calculation, and vice versa, is not found 
to describe the data well, enforcing the idea that the en- 
hancement is substantially and fundamentally different 
for these two elements. Since V and Cr are electroni- 
cally very similar, exhibiting the same body-centered cu- 
bic structure, the largest difference between the two is in 
their band filling and FS. In V, Ep lies close to a peak 
in the d density of states with appreciable {^ 20 %) p- 
character, leading to a total number of states at £^f of 
N{E-p) = 23.80 states per Ry per atom. In paramagnetic 
Cr, on the other hand, the additional electron places Ep 
in a valley between the bonding and anti-bonding d-states 
with N{Ey) = 9.52 states per Ry per atom. It follows 
that the number of electrons that are capable of screening 
the positron impurity (and thus lead to the enhancement 



of the annihilation rate) in V and Cr is very different, 
and cannot be captured by considerations of the elec- 
tron density or energy alone. However, such a concept 
does provide a route to understanding the different SD 
enhancement models in V and Cr, and the different con- 



stants b/a and c/a in Eq. 10 that are required to explain 
the data. 

In order to test such a correction to the enhancement, 
we apply a MS-type enhancement to the SBA model, 
which is then fitted to the experimental data. According 
to Kahana's theory, the quadratic part of the enhance- 
ment parameterization is fairly weak, with c/a ~ 0.138 
for metallic densities, and can be well approximated by 
just a linear component {b/a). Here, we adopt just this 
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linear energy enhancement, and set c/a = in Eq. 
leaving just a single fitting parameter that is capable of 
adjusting the shape of the computed RMD: 



7SBA-MS = 7sba[1 + {b/a)e\. 



(11) 



The results of such a model, which we refer to as SBA- 
MS enhancement, are found to enormously improve the 
agreement between data and theory for all materials 
(see Table VII), leading to x?ed parameters that ap- 



proach the SD model investigated in Section HI More- 
over, for the three transition metal elements addressed 
in this manuscript, this linear component of energy en- 
hancement is found to sc ale w ith the density of states 
at the Fermi level (Table | VII ) , providing an empirical 
model for the enhancement of d-band elements and com- 
pounds. For s-p electron metals, the enhancement is 
found to more closely resemble Kahana's parameters, and 
two regimes emerge - Kahana's prediction for s-p simple 
metals, and a A^(^F)-dependent set of parameters for d 
electron metals. It is interesting to compare our results 
for an MS-type enhancement to those applied by Mat- 
sumoto and Wakoh for Cr [21 , in which they obtain a 
factor b/a ~ 0.15 in a non-optimized approach, very close 
to our 0.20. On the other hand, Genoud [22 , employing 
an MS-type enhancement independently for 5, p and d 
electrons in V, obtained b/a = 0.1 — 0.2 for s and p elec- 
trons, also in a non-optimized way, which is somewhat 
smaller than the optimum b/a = 0.70 that we find. 

Armed with such a model for the enhancement, we 
can now assess its validity for another metal, specifically 
the d electron metal Mo, for which we have just a single 
[110] projection available, insufficient to permit a reliable 
fit of the ji parameters of the SD model. Instead, we ap- 
ply the SBA-MS model to the data, in which the 7/'s are 
computed from the partial annihilation rates of state / 
and the energy-dependent enhancement is provided from 
N{Ef) = 7.6 states / Ry / atom and our preceding fit. 
For comparison, we also compute the IPM and JS RMD. 
Both the JS and SBA models of enhancement, by them- 
selves, offer negligible improvement over the IPM, which 
already provides a reasonable description of the data, and 
the application of the SBA-MS model improves the agree- 
ment only modestly by ~ 4 %. However, lending freedom 
to the linear component b/a is not found to provide any 
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b) 



^iwu 




c) 



d) 




FIG. 8: (color online) Comparison between experimental data 
and the corresponding theoretical quantities in the raw IPM 
and rigid-band fit SD models of enhancement for (a) Cr and 
(b) Mo. The RMD is shown along a path in the BZ in (c) 
and (d) for Cr and Mo respectively. 



additional improvement, emphasizing that the original 
IPM calculation was already satisfactory. 

One of the unresolved questions of 2D- AGAR in tran- 
sition metals is why the k-space density of Gr and Mo 
(projected along the [110] direction) appear so different, 
despite the apparent similarity of their isoelectronic and 
isostructural FS topology [4T (see Fig.|8|i,b). In Ref. |44j . 
maximum-entropy filtering techniques were employed to 
assess the FS breaks in both distributions, ruling out the 
FS topology as an explanation; the question of whether 
positron effects or consequences of the proximity to mag- 
netic structure in Gr are to blame were left open and 
have remained so despite several efforts to resolve the is- 
sue, both experimentally and theoretically [53H55] . Here, 
we are able to solve this issue, which stems from a strong 
over-estimation of the enhancement near FN in Gr, pre- 
viously highlighted in Fig. (S) In Fig. [8|3,d, the RMD 
along a path in the BZ is shown for both Gr and Mo for 
the 2D- AGAR data, and the IPM and SBA-MS models 
of enhancement. It is clear that the IPM (which resem- 
bles the JS model) looks similar for both metals, elimi- 



nating positron wavefunction effects (which are included 
in the IPM) as responsible for the strong difference in 
the data. On the other hand, the SBA-MS prediction 
(which closely resembles our measured SD model), ac- 
counts for the data very closely, unambiguously estab- 
lishing enhancement effects as the key. 



VI. CONCLUSIONS 

We have presented a detailed investigation of the 
positron enhancement factor for several metals, provid- 
ing a quantitative measurement of the state-dependence 
of the enhancement. By combining this with a rigid 
shift of the energy bands, we demonstrate that, when 
the band structure is optimized to 2D-AGAR measure- 
ments, the precise location of the Fermi breaks in k- 
space are sensitively dependent on the accuracy of the 
enhancement model used in the calculation. Further- 
more, we show that, by employing a state-dependent 
model for the enhancement, much improved agreement 
between the 'tuned' calculation and high-precision quan- 
tum oscillatory data can be obtained. In particular, for 
Gr our positron measurements yield a nesting vector that 
is in excellent agreement with neutron measurements of 
the spin-density wave ordering vector, with an estimated 
accuracy better than 0.5% of the BZ. Although alloys 
have not been investigated here, this approach also al- 
lows for the contribution from different atomic sites to 
be separated in the experimental data, allowing a deter- 
mination of the fraction of annihilations from each indi- 
vidual element's hybridized wavefunctions (for example, 
see Ref. [56]). 

Comparisons of our (empirical) model with other pop- 
ular models of the enhancement have been made, partic- 
ularly with the ab initio state-dependent model of Bar- 
biellini, Alatalo and co-workers [5 , for which a semi- 
empirical energy-dependent correction is proposed that 
is found to bring the theory into much better agreement 
with the data. Such a combined model therefore pro- 
vides an accurate model for the enhancement in momen- 
tum density measurements, such as those of 2D-AGAR 
or coincidence Doppler broadening techniques. 
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